Two-component model of a spin-polarized transport 
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Effect of the spin-involved interaction of electrons with 
impurity atoms or defects to the transport properties of a two- 
dimensional electron gas is described by using a simplifying 
two-component model. Components representing spin-up and 
spin-down states are supposed to be coupled at a discrete set 
of points within a conducting channel. The used limit of the 
short-range interaction allows to solve the relevant scattering 
problem exactly. By varying the model parameters different 
transport regimes of two-terminal devices with ferromagnetic 
contacts can be described. In a quasi-ballistic regime the 
resulting difference between conductances for the parallel and 
antiparallel orientation of the contact magnetization changes 
its sign as a function of the length of the conducting channel 
if appropriate model parameters are chosen. This effect is in 
agreement with recent experimental observation. 



Spin-polarized transport in two-dimensional electron 
systems has been a field of growing interest during the 
last years. Typically, the experiments are performed 
using a two-terminal device with ferromagnetic metal 
contacts. A spin-polarization of the injected current is 
expected from the different densities of states for spin- 
up and spin-down electrons in the ferromagnetic source. 
This leads to a spin dependent interface-resistance, which 
also exists at the interface of the second ferromagnetic 
contact, the drain. Together with spin-involved scatter- 
ing processes in the studied electron system this should 
result in a conductance which depends on the relative 
magnetization of the two contacts jlj . 

The quantum mechanical nature of spin places it out 
of reach of many of the forces in a solid and the ori- 
entation of a carrier's spin can be very long-lived. The 
conductance G'' of a two-terminal device with parallel 
orientation of magnetic moments of the contacts is thus 
expected to be higher than the conductance G**- for the 
case of antiparallel moment orientation However, 
just the opposite results have been reported recently 0] 
for a two-dimensional electron gas confined in an InAs 
channel with the permalloy source and drain. It has been 
found that an ensemble average of the conductance differ- 
ence G^T — G?tJ. decreases as function of the channel length 
reaching negative values in a quasi-ballistic regime when 
the electron mean free path l e becomes comparable with 
the channel length. 

In the absence of magnetic impurities the natural can- 
didate for spin dephasing and precession effects is spin- 
orbit coupling to impurity atoms or defects. General the- 
oretical approach to its description is to use the contri- 
bution to the Hamiltonian which stems directly from the 



quadratic in v jc expansion of the Dirac equation 



Hso 



T VV(r)-(crxp) , 



(1) 



where a = {a x , a y , a z } denotes Pauli matrices, V(r) is a 
background potential, V stands for the spatial gradient 
and m is the electron mass. 

Influence of the spin-orbit interaction on the electron 
transport properties of two-dimensional mesoscopic sys- 
tems has been studied since the early 1980's. At that 
time it was found that it is responsible for so-called an- 
tilocalization effect Q|. Later the attention has been 
turned to the effects caused by a Rashba term J5,|| in 
two-dimensional ]^[Tc|] and quasi-one-dimensional sys- 
tems pl|-[l5|. Realistic transport theory for fully quan- 
tum coherent systems including the spin-orbit coupling to 
the impurities or defects has not yet been reported. The 
problem becomes complicated even if electron motion is 
restricted to the two-dimensional space. In general, the 
spin-orbit term, Eq.(Q), turns the problem back to three 
dimensions. 

The goal of this paper is to reveal those features of the 
transport properties which can be caused by the spin- 
orbit interaction induced by a scattering potential. In the 
interesting case of a quasi-ballistic regime, which shows 
chaotic features, it is very difficult to estimate deviation 
from the exact solution caused by any used approxima- 
tion. For this reason we have employed simplifying two- 
component model with point interaction for which exact 
solution, including fully the quantum coherence, can be 
found. Non-zero spin-orbit coupling is assumed to be as- 
sociated with short-range scattering potentials only. Al- 
though the treatment is far from a realistic transport 
theory, it might be useful to understand some mysteries 
of the recent experimental observation. 

Free electron system is a typical two-component sys- 
tem if the electron spin is taken into account. If there are 
no spin-involved forces, electron states are represented by 
plane waves exp(ikr) with k being a wave vector. Ori- 
entation of the electron spin is given by the quantum 
number s z = ±a and the electron system can be splittcd 
into two independent subsystems, each of them composed 
of electrons having the same spin orientation. However, 
any perturbation of the background potential can cause a 
coupling between subsystems due to non-zero spin-orbit 
term defined by Eq.(|l]). 

Let us first consider a single scattering potential acting 
on a two-dimensional electron gas within a finite region 
of a radius tq. In accord with standard scattering theory 
an incoming wave belonging to one particular subsystem, 



say of the spin up states, can be scattered into states be- 
longing to both subsystems. In the short-range limit, 
krQ << 1, only s-part of the incoming wave gives non- 
zero contribution to the scattering process. For given en- 
ergy E = h k 2 /2m, (k = |fc|), the corresponding solution 
of the radial Schrodingcr equation has two components, 
^l(r) and \&i(r). Outside the scattering region they can 
be written as follows ]l6| , [l7t 



* T (r) = Jo(fcr) + a(k)H^\kr) 
#i(r) = b(k)H^(kr), 



(2) 
(3) 



where Jo(z) denotes the cylindrical Bessel function and 

Hankel functions Hq (z) represent scattered outgoing 
waves that for large arguments have the following asymp- 
totic form 



H^\kr ->oo) - ^-== 



1 
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(4) 



Taking into account the time reversal symmetry and as- 
suming that the system is invariant with respect of the 
subsystem interchange, the amplitudes a(k) and b(k) in 
the short-range limit have to be of the following general 
form Q 



a(k) 



[1 
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b(k) = -C 
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a{k) , 



(5) 
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where A and C are real model parameters. If C is chosen 
to be zero, the parameter A represents a strength of the 
scattering process within one particular subsystem and 
for a potential well of the radius ro it takes the value 
A = In ro. Non-zero values of the parameter C give rise 
to a spin- flip process. 

There arc two relevant physical quantities characteriz- 
ing scattering event: the total scattering cross-section 



(To = a(k) 2 + b(k) 2 

and the spin-flip probability t^-l 

b(k) 2 
a{k) 2 + b(k) 2 



(7) 



(8) 



Note, that the assumption of the system invariance with 
respect of the subsystem interchange leads to the inde- 
pendence of ao and t^-l on the spin orientation of the 
incoming electron. This assumption has been used for 
the sake of simplicity despite of the fact that it need not 
to be satisfied in real systems, e.g. due to a Rashba term. 

The scattering problem for a two-dimensional strip 
with a finite number of short-range scatterers, as sketched 
in Fig. 1, can be solved exactly. The detail analysis in 
the case of a one-component system has already been 
reported |18| and generalization to a two-component sys- 
tem is straightforward. For simplicity we have assumed 



that all scatterers are identical, i.e. they give the same 
scattering cross-section ao and spin-flip probability t'^ if 
they would be placed alone within the two-dimensional 
space. Scattering matrix has been obtained numerically 
for a given configuration of point scatterers randomly 
distributed within a strip region of the length L. 
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FIG.l. Scheme of the scattering process in a two-component 
system. Upper and lower strips of the width w represent the spin- 
subsystems. Scatterers, black points, serve also as connection 
points between subsystems giving rise to a spin-flip processes. 
Thick full and dashed lines represent an incoming and outgoing 
waves, respectively. 

Spin-dependent transport properties are determined 
by partial transmission coefficients representing transi- 
tion between left and right subsystems of asymptotic 
spin-up or spin-down states. They are defined as the 
sum of transmission probabilities over all relevant modes 
of asymptotic states. To simplify the description by ex- 
cluding the quantum fluctuations from our consideration 
we have used configurationally averaged values of the par- 
tial scattering coefficients to define 2x2 transmission and 
reflection matrices T and R, respectively: 



R 



i?TT R n 



(9) 



For the considered symmetrical system T^T = j>U ) jt.1 = 
= and R n = R 1 ^ . 




FIG. 2. The total transmission coefficient T as function of the 
scattering region length L for several values of the spin-flip prob- 
ability: t n = 8.1 X 10" 3 (full line), $Ti = 4.2 X 10" 3 (dashed- 



dotted line), t n 
(dashed line) . 



1.1 X 1(T 3 (dotted line), *tl = 0.0005 X 10" 
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In Fig. 2 and Fig. 3 the dependence of transmission 
coefficients on the length L of the scattering region is 
shown for different spin- flip probabilities V^. The used 
energy corresponds to 31 occupied subbands. Concen- 
tration of scatterers (750/w 2 ) and the scattering cross- 
section (To = 0.1217 were held fixed. Probability of an 
injected electron to cross scattering region increases with 
increasing spin-flip probability as shown in Fig. 2 where 
the total transmission coefficients T = 2(T^ + T^i) are 
plotted. This tendency is in agreement with expected 
antilocalization effect M . 
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FIG. 3. The partial transmission coefficients T^T (full line) and 
T^i (dashed line) as function of the scattering region length L 
for several values of the spin-flip probability: a) = 8.1 X I0 — 3 ; 
b) tt-t = 0.2 x f(T 3 ; c) ttJ. = 0.0005 x f0~ 3 . 

The more interesting is the dependence of partial 
transmission coefficients. For some values of the spin-flip 
probability and lengths L, becomes less than T'*, 
as can be seen in Fig. 3. It means that the polarization 
of the transmitted current has opposite orientation than 
the polarization of the injected current. This surprising 
result we ascribe to the non-trivial quantum coherence 
in two-component systems. Wave interference leading to 
the weak localization is one-component effect. It take 
place within one particular subsystem only. In the con- 
sidered case of a weak spin- flip process (t^ << 1) it 
becomes dominant within the subsystem with incoming 
waves while the localization effect within second subsys- 



tem do not need to be for given length L still well devel- 
oped. 

The above described effect disappears if the scattering 
cross-section is substantially enlarged. The localization 
becomes dominant effect of the scattering process in the 
both subsystems and inequality T" > T*^ remains valid 
for all lengths of the scattering region. 

The device conductance of a two-component quantum 
system is determined by a matrix of partial transmission 
coefficients, Td GV , 

G=£(l,l)Tdev(j) ; T dev = ^ (1Q) 

which depends on the properties of ferromagnetic con- 
tacts and their interfaces with the two-dimensional elec- 
tron gas. To estimate their effect we have used the idea of 
polarization filters [l4| . The source and the drain are con- 
sidered to be standard reservoirs and all spin-dependent 
effects are modeled by filters placed within the asymp- 
totic region of ideal leads. If the coherence is supposed 
to be completely destroyed at the filter boundaries the 
conductance can be expressed as a function of the al- 
ready defined coefficients T'' , R*' and describ- 
ing scattering process of the same device without filters. 

Experiments on spin-injection into a two-dimensional 
systems usually show a large interface resistance. For 
this reason non-zero probabilities, a and 0, of spin-up 
and spin-down electrons to be reflected by the filter will 
be considered. For the sake of the simplicity we assume 
that reflected electrons will be equally distributed be- 
tween available quantum channels without any change of 
their spin orientation. In this case the filtering effect can 
be described by 2 x 2 diagonal matrix 

where indices s and d represent the source-filter and the 
drain-filter, respectively. 

For the case of N available quantum channels (sub- 
bands) within each subsystem, the above described 
model leads to the following expression for the transmis- 
sion matrix Td ov entering the Eq. (|To|) 

T dev = (1 - F d ) NM d T K rf , s NM S (1 - P.) (12) 

where N stands for the product of N and unit matrix 1. 
The effect of multiple reflections between filters and the 
scattering region is represented by matrices M s and 

M s = (N - FsR)- 1 ; M d = (N - RF^ 1 , (13) 

and 

K d . s = [1 - MsF.TFrfMrfT]- 1 . (14) 

Device conductance, Eq.(|l0|) depends on reflection 
probabilities a Std and (3 S:d modelling the effect of ferro- 
magnetic contacts. For the case of the parallel orientation 
of the contact magnetization the conductance G'' can be 



obtained by setting a s = ad and (3 S = (id- To get G^l 
for the antiparallcl contact magnetization a s = [3d and 
(3 S = ad have to be used. While the conductance strongly 
depends on the used values of reflection probabilities, the 
sign of the conductance difference AG = — G^, is 
not affected. 
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FIG. 4. Relative conductance change AG/Go as function of the 
scattering region length L for several values of the spin-flip prob- 
ability: *ti = 8.1 X (full line), *Ti = 1.1 x 10" 3 (dashed- 
dotted line), tT-l = 0.2 X 10~ 3 (dotted line), *ti = 0.0005 X 10" 3 
(dashed line). In the inset crosses represent experimental data 
obtained by Hu et al. and full line is the result of the model 
calculation for the following parameters: N = 173, a s = 0.01, 
/3 S = 0, scatterer concentration 1500/w 2 , <ro = 0.1341 and 
tTi = 4.8 x io~ 3 . 



In Fig. 4 the relative conductance change 

AG _ G TT - G n 
G7 = GTT + GTT 



(15) 



as a function of the scattering-region length L is shown 
for the case of ideal filters, a s = and (3 S = 1. It corre- 
sponds to injection of fully polarized current and vanish- 
ing interface resistance. All other used model parameters 
are the same as that for transmission coefficients pre- 
sented in Fig. 2 and Fig. 3. The model parameters of the 
scatterers, <7o and t^, and their concentration have been 
chosen to give pronounced minimum in the dependence 
AG/Go on L. 

To model the device studied by C. M. Hu et al. with 
hundreds occupied subbands it is necessary to perform 
calculation for much larger energy and take into account 
boundary resistance between ferromagnetic contacts and 
two-dimensional electron gas by using non-zero value for 
a s . Note that a s ~ 0.01 lowers the values of the relative 
conductance change approximately hundred times and 
nearly quantitative agreement with the measured data 
can be reached as shown in the inset of the Fig. 4. 

The presented two-component model does not allow to 
consider asymmetry of scattering processes which leads 
to the spin-Hall effect [Q. However, if the polarization 
vectors of the contact magnetization lies within the plane 
of the electron gas and they are perpendicular to the 



applied current, as in the case studied by C. M. Hu et al. 
||, this effect is suppressed. 

The main result of the described model is that in 
mesoscopic disordered systems the quantum coherence 
affected by spin-flip processes can lead to the higher con- 
ductance of two-terminal devices with antiparallel con- 
tact magnetization than that for parallel configuration 
as observed in recent experiments by by C. M. Hu et 
al. ||. Similar effects can also be expected for other two- 
component systems. For example in double-layer systems 
the defects of the barrier separating two-dimensional elec- 
tron gases can act as their connection points. Instead of 
a spin-flip process there will be a real electron transition 
between electron layers and in particular cases the con- 
ductance measured between contacts at different layers 
could be larger than that between contacts of the same 
electron layer. 
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